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ABSTRACT 

The cooling-flow problem is a long-standing puzzle that has received considerable recent 
attention, in part because the mechanism that quenches cooling flows in galaxy clusters is likely 
to be the same mechanism that sharply truncates the high end of the galaxy luminosity function. 
Most of the recent models for halting cooling in clusters have focused on AGN heating, but 
the actual heating mechanism has remained mysterious. Here we present a framework for AGN 
heating derived from a Chandra survey of gas entropy profiles within cluster cores. This set of 
observations strongly suggests that the inner parts of cluster cores are shock- heated every ~ 10 8 
years by intermittent AGN outbursts, driven by a kinetic power output of ~ 10 45 ergs -1 and 
lasting at least 10 7 years. Beyond ~ 30 kpc these shocks decay to sound waves, releasing buoyant 
bubbles that heat the core's outer parts. Between heating episodes, cooling causes the core 
to relax toward an asymptotic pure-cooling profile. The density distribution in this asymptotic 
profile is sufficiently peaked that the AGN shock does not cause a core entropy inversion, allowing 
the cluster core to retain a strong iron abundance gradient, as observed. 



Subject headings: galaxies: 
galaxies: clusters 

Introduction 



clusters: general — galaxies: evolution — intergalactic medium — X-rays: 



The gaseous cores of galaxy clusters have de- 
fied interpretation for almost three decades. Ini- 
tially, gas in the cores of clusters was thought to 
cool, condense, and flow toward the center, as long 
as the time required for the central gas to radi- 
ate away its thermal energy was less than the age 
of the universe (Fabian & Nulsen 1977; Cowie & 
Binney 1977; Mathews & Bregman 1978). Most 
galaxy clusters satisfy this condition, and those 
that do have traditionally been called cooling-flow 
clusters. Early estimates of their mass accretion 
rates ranged as high as 10 2 — 10 3 M & yr (Fabian 
1994). The problem with this interpretation of the 
data was that the mass sink for all this supposedly 
cooling and condensing gas has never been found 
(Braine & Dupraz 1994; McNamara & Jaffe 1994; 
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O'Dea et al. 1994, 1998; Voit & Donahue 1995; 
Donahue & Voit 2004). 

Spectroscopic X-ray observations have now es- 
tablished that the rates of cooling and condensa- 
tion in galaxy clusters are much lower than those 
originally inferred from the imaging observations 
(Peterson et al. 2001, 2003; Tamura et al. 2001). 
There is little evidence for emission lines from gas 
cooling much below about half the ambient tem- 
perature of the hot gas filling the cluster. The 
cooling process seems to have stopped in its tracks, 
implying that a compensating heat source resup- 
plies the energy radiated by the gas in cluster 
cores. 

Supermassive black holes at the centers of 
galaxy clusters are an attractive candidate for 
supplying that heat energy (e.g., Binney & Tabor 
1995). Clusters suspected of harboring cooling 
flows are statistically more likely to contain radio 
emission indicative of recent black-hole activity 
(e.g., Burns 1990), and that radio emission of- 
ten coincides with cavities in the X-ray emitting 
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plasma, clearly showing that outflows from the 
nuclei of galaxies are interacting with the core gas 
in clusters (e.g., McNamara et al. 2000; Fabian 
et al. 2000). Measuring the work needed to in- 
flate these cavities reveals that each outburst of 
energy from an active galactic nucleus (AGN) in- 
troduces ~ 10 58 — 10 61 erg into the intracluster 
medium (Birzan et al. 2004; McNamara et al. 
2005), which would be sufficient to stop the cool- 
ing if such outbursts occurred every ~ 10 8 years 
or so, supplying a time-averaged energy input rate 
- 10 43 - 10 45 erg s- 1 (Churazov et al. 2002; David 
et al. 2001). 

While AGN heating is certainly plausible, it 
has not been clear how the outburst energy gets 
transfered to the intracluster gas (cf. Babul et al. 
2002). Suggestions that the AGN outflow heats 
the ambient gas by driving a powerful shock (e.g., 
Soker et al. 2001) were dismissed when early Chan- 
dra observations showed cool rims rather than hot 
ones around the X-ray cavities, motivating inves- 
tigations into gentler heating mechanisms (Chura- 
zov et al. 2000, 2001; Begelman 2001; Ruszkowski 
& Begelman 2002; Quilis et al. 2001; Reynolds 
et al. 2002; Briiggen & Kaiser 2002; Dalla Vecchia 
et al. 2004). The cool rims now seem to consist 
of low-entropy gas transported to larger radii by 
the buoyantly rising cavities, but there are other 
reasons why gentle heating mechanisms have been 
preferred. One is that clusters with short central 
cooling times have iron abundances that strongly 
decrease with radius (e.g. De Grandi et al. 2004). 
An impulsive injection of energy at the center of 
the cluster would appear to destroy such a gradi- 
ent by causing the cluster core to convect, thereby 
mixing the iron more uniformly (but see Omnia 
et al. 2004). Another is that the abrupt jumps in 
gas pressure expected of shocks were not obvious 
in much of the early Chandra data; however, de- 
tecting a localized region of hot gas can be quite 
difficult if emission from lower-temperature gas is 
present along the same line of sight (e.g., Mazzotta 
et al. 2004). 

Some recent X-ray observations are suggesting 
that the shock-heating hypothesis may have been 
prematurely discarded. Deep Chandra images of 
Hydra A (Nulsen et al. 2004) and MS0735+7421 
(McNamara et al. 2005) have revealed that AGN 
activity can drive shock fronts to distances ~ 200— 
300 kpc from the center, indicating that a kinetic 



power output ~ 10 45 — 10 46 ergs _1 has persisted 
for ~ 10 8 years in those clusters. Cases like these 
may be rare, but objects that currently show ra- 
dio emission indeed seem systematically different 
from ones that do not. Croston et al. (2004) have 
reported that groups of galaxies with radio emis- 
sion have systematically lower X-ray luminosities 
and higher X-ray temperatures than those with- 
out radio emission. They interpret this difference 
as transient spikes in temperature because the ra- 
tio of optical to X-ray luminosity shows no such 
systematic difference. 

Another clue comes from Chandra observa- 
tions of two radio-quiet clusters by Donahue et al. 
(2005b). Both clusters were classified as strong 
cooling flows by (Peres et al. 1998), but their cen- 
tral cooling times (~ 10 9 years) are much longer 
than the central cooling times of cooling-flow clus- 
ters with current radio activity (~ 10 8 years). 
Either the gas in these clusters has never cooled 
and condensed, obviating the need for feedback, or 
they experienced such a dramatic feedback episode 
sometime in the past that they have not required 
any additional feedback for ~ 10 9 years. The fact 
that signs of AGN feedback appear only in those 
cooling-flow clusters that currently require feed- 
back strengthens the case for AGN heating as the 
mechanism that limits cooling. 

Here we revisit shock heating in the context 
of a Chandra sample of classic cooling-flow clus- 
ters (Donahue et al. 2005a). Profiles of gas en- 
tropy in the cores of these clusters show a striking 
consistency that stands as a strong constraint on 
any model purporting to explain the cooling-flow 
problem. Section 2 assesses the clues these pro- 
files hold about the cooling and feedback mech- 
anisms in cluster cores, while Section 3 outlines 
a framework for AGN heat input consistent with 
those clues. The picture we arrive at from an- 
alyzing these entropy profiles agrees in many re- 
spects with the AGN-heating models inspired by 
the X-ray cavities and is aligned with the sugges- 
tion of Binncy (2005, and references within) that 
episodic outflows from AGNs quench cooling and 
condensation in clusters, limit the maximum lu- 
minosity of galaxies, and regulate the growth of 
black holes at their centers. We show that out- 
bursts of ~ 10 45 ergs -1 occurring on a ~ 10 8 year 
timescale are sufficient to maintain the observed 
core entropy profiles. However, unusually pow- 
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crful outbursts like those seen in Hydra A and 
MS0435+7241, which are rarer and longer-lasting 
than those needed to maintain the observed pro- 
files, may also play a role in extending the time 
between outbursts and in elevating intracluster en- 
tropy beyond the cores of clusters. Section 4 sum- 
marizes our results. 

2. Core Entropy Profiles 

We focus here on the entropy content of intr- 
acluster gas rather than on temperature and gas 
density individually because entropy offers more 
direct insights into the processes that add and re- 
move thermal energy from the gas. Adding heat 
energy to a cluster does not simply raise the gas 
temperature. Instead, the heat input causes the 
cluster's gas to expand within its potential well, 
and the temperature it achieves after relaxing to 
equilibrium depends primarily on its entropy dis- 
tribution and the shape of that potential well (e.g., 
Voit et al. 2002; Voit 2005). In this paper we 
will quantify entropy with the adiabatic constant 
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K = kTn e , because this quantity provides the 
most direct link between X-ray observations and 
entropy content. 

Figure 1 shows nine core entropy profiles of 
cooling-flow clusters taken by Donahue et al. 
(2005a) from the Chandra archive: 2A0335+096, 
Abell 133, Abell 262, Abcll 496, Abell 1795, 
Abell 2029, Abell 2052, Hydra A, and PKS 
0745+191. The data were azimuthally averaged 
in annuli centered on the cluster core and then de- 
projected, giving n e and T as functions of radius, 
which were then combined into K(r) profiles. The 
similarity among the profiles immediately suggests 
that some sort of quasi-steady process stabilizes 
the intracluster gas. There are no large entropy 
inversions outside of 10 kpc, and the apparent 
entropy inversion seen at ~ 9 kpc in one of the 
clusters (Abcll 2052) arises primarily from non- 
axisymmetric structure owing to its X-ray cavi- 
ties. 

Also notable is the fact that entropy at a given 
physical radius is quite similar from cluster to clus- 
ter, even though the clusters range in tempera- 
ture from 2.2 keV to 7.4 keV. That form of sim- 
ilarity can be understood if the entropy profiles 
of clusters scale with the characteristic halo tem- 
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perature T$ according to K oc TV , as suggested 
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Fig. 1. — Core entropy profiles of nine classic cooling-flow 
clusters from Donahue et al. (2005a). The dotted lines show 
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the gas entropy (K = kTn e ) profiles derived from the 
azimuthally averaged and deprojected density and temper- 
ature profiles of these clusters as a function of radius. The 
profiles look strikingly similar considering that the cluster 
temperatures in this sample range from 2.2 keV to 7.4 keV. 
Similarity in the entropy levels of the outer regions is con- 
sistent with the T 2 ' 3 scaling of the outer entropy profiles 
of clusters seen by Ponman et al. (2003). At larger radii 
the profiles seem to asymptotically approach the solid line 
showing the expected entropy profile for a 5 keV cluster 
model from Voit et al. (2002) in which radiative cooling 
acts without triggering feedback. Simply adding a constant 
10 keV cm 2 to the asymptotic model gives the thick dashed 
line, which is a good representation of the observed profiles 
within about 50 kpc (assuming Hq = 70kms _1 Mpc -1 ). 

by the data of Ponman et al. (2003). Donahue 
et al. (2005a) find that the dependence of entropy 
on radius in these clusters is K{r) oc r a with 
a « 0.9 — 1.3. Assuming that a = 1.1 and that 
the virial radius of a cluster's halo is oc T^ 2 , one 

finds that K (r) oc T 2/3 (r/T^ 2 ) 1 - 1 oc T°- 12 r 1A . In 
other words, the K(r) profiles depend only weakly 
on temperature, as long as the entropy at a given 
fraction of the virial radius is oc T 2/3 . It is not yet 
clear whether this altered entropy scaling stems 
mostly from cooling and supernova feedback dur- 
ing the epoch of galaxy formation or whether it 
requires additional AGN heating at lower redshifts 
(e.g., Bryan & Voit 2005) 

The thick solid line in Figure 1 shows the core 
entropy profile one would expect for a 5 keV clus- 
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ter in which radiative cooling operates without any 
feedback, according to Voit et al. (2002). Their 
model started with the baseline entropy profile 
expected from gravitational structure formation 
alone and then allowed that gas to cool for a Hub- 
ble time within a Navarro-Frcnk- White potential 
(Navarro ct al. 1997). The inner temperature pro- 
files predicted by that pure cooling model (the Kb. 
model in the lingo of Voit et al. (2002)) were shown 
to be quite similar to those of clusters suspected to 
harbor strong cooling flows, suggesting that cool- 
ing indeed plays an important role in determining 
their core entropy structure (see also McCarthy 
et al. 2004). 

Donahue et al. (2005a) show that cooling-flow 
clusters do not actually relax to the pure-cooling 
profile. Instead, they appear to asymptotically 
approach it at ~ 100 kpc. Meanwhile, the in- 
ner entropy profiles of these clusters flatten out 
at the ~ lOkeVcm 2 level, a clue to the form of 
heat input that hinders cooling and condensation. 
Assuming pure free-free cooling, one can write a 
simple relation between entropy cooling time, and 
gas temperature: 

The central flattening of these clusters' entropy 
profiles therefore suggests that they will not re- 
quire feedback to halt cooling until ~ 10 8 years 
from now. In that respect, the clusters in the 
Donahue et al. (2005a) sample follow a trend first 
noticed in Hydra A by David et al. (2001). Heat- 
ing would therefore seem to be episodic (Kaiser 
& Binney 2003). However, in order to maintain 
the quasi-steady nature of these profiles, the max- 
imum entropy input in a typical heating event 
cannot be much greater than 10 — 20keVcm 2 . 
Otherwise the cluster core would linger at much 
higher entropy levels, because it takes much longer 
to cool from a high-entropy state, owing to the 
steeper-than-linear dependence of t c on K. The 
dashed line in Figure 1 shows that simply adding 
a constant 10 keV cm 2 entropy pedestal to the pure 
cooling model produces a reasonable facsimile to 
the observed profiles out to > 30 kpc, another clue 
to the mode of heat input. 

The structure of these cluster cores from 
<~ 10 — 100 kpc can be approximated with a sim- 
ple power-law model that will be useful when 



we consider heat input. Figure 1 shows that 
the entropy profiles in this range are approxi- 
mately K(r) w 150keVcm 2 (r/lOOkpc), and the 
temperature gradients at these radii in cooling- 
flow clusters typically have a power-law slope 
d\iiT/d\nr w 1/3 (e.g., Voigt & Fabian 2004; 
Donahue et al. 2005a). Combining these facts, 
we can reverse-engineer the observed density pro- 
file n e (r) w 6.1 x 10~ 3 cm" 3 T^ Q0 (r/lOOkpc)" 1 , 
where T^ioo is the gas temperature 100 kpc from 
the cluster center in units of 5 keV. The product 
rp(r) « 3.7 x 10~ 3 gcm~ 3 (q is therefore ap- 
proximately constant in this radial range, and is 
consistent with the observations of Donahue et al. 
(2005a). 

Before we discuss heat input, we have a com- 
ment on this power-law model. Notice that the 
gas temperature drops by a factor of two from 100 
kpc to 10 kpc. Inside of 10 kpc, the gas tem- 
perature must then level off, and may even rise a 
little, if the gas within 10 kpc is isentropic. Thus, 
a heating mechanism that maintains an inner en- 
tropy plateau at ~ lOkeVcm 2 for long periods of 
time would account for the spectroscopic finding 
that gas temperatures in the centers of cooling- 
flow clusters drop to about half the ambient tem- 
perature, even in a quasi-static configuration in 
which there is little net cooling. 

3. A Framework for Heat Input 

Consider how an energetic outburst from an 
AGN heats gas with this power-law configura- 
tion. An outflow with a kinetic power output 
-Pkin ~ 10 45 ergs -1 shock-heats the innermost por- 
tions of the core. Even if the flow is bipolar, an el- 
lipsoidal cocoon of shocked intracluster gas devel- 
ops around the outflow, distributing heat through 
47r steradians (e.g., Scheuer 1974; Begelman & 
Cioffi 1989). If the outflow shuts down after ~ 10 7 
years, the energy introduced, amounting to a few 
times 10 59 erg, continues to drive the shock until 
it becomes subsonic. The hot bubbles inflated by 
the outflow then become buoyant, floating away 
from the center and gradually thermalizing the 
PdV work they exert on their surroundings (e.g., 
Churazov et al. 2001; Begelman 2001). Here we 
show that the entropy boosts imparted by a series 
of such episodic heating events are broadly consis- 
tent with the observed entropy profiles of cooling- 
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flow clusters. 

3.1. Outflow-Driven Shock Heating 

We can estimate the entropy introduced dur- 
ing the outflow-driven stage using the scaling be- 
havior of an outflow-driven bubble and an entropy 
jump condition. The shock velocity of such a bub- 
ble scales as v w fp(Pkin/ pr 2 ) 1 ^ 3 , where fp is a 
structure factor of order unity that depends on the 
outflow geometry and the preshock density profile 
(e.g., Ostriker & McKee 1988). In order to assess 
the entropy created by this shock, we apply the 
approximate entropy jump condition 



AK = Ko — Ki 



pm p v 
3(4n e ) 2 / 3 



-0.16 If i 



(2) 



where K x and K 2 are the pre- and post-shock en- 
tropy levels, respectively, and n e is the preshock 
electron density (Voit et al. 2003). For the pur- 
poses of this estimate, we will ignore the insignif- 
icant 0.16 ifi term. The entropy increase AKp 
induced by such an outflow-driven shock in a 
medium with a constant value of rp(r) is then in- 
dependent of radius, with AKp oc Pk!n( r P)~ 4 ^ 3 ■ 
Its physical value in the power-law core model of 
§2 is 



AK P « 23 keV cm 2 P 2 / 3 f P 



(3) 



where P 45 = P kin /10 45 erg s" 1 . 

This mode of entropy injection has the very at- 
tractive feature of introducing a constant entropy 
boost throughout the shock-heated zone, as sug- 
gested by the observations. Outflow-driven shocks 
create no entropy inversions and therefore do not 
induce the kinds of mixing that would destroy the 
iron gradients that have built up in cooling-flow 
clusters. After the shock passes and the gas set- 
tles back into equilibrium, the stratification of the 
cluster core remains the same because the order- 
ing of specific entropy among gas parcels has not 
changed. If this is indeed the primary mode of en- 
tropy injection within the centers of clusters, then 
the observed ~ lOkeVcm 2 entropy plateau im- 
plies a kinetic power output ~ 10 45 ergs -1 during 
each heating event. Such an outflow-driven shock 
heats a region extending to a radius 



rp « 16kpc 



AK f 



lOkeVcm 2 



3/8 



r 3/8 ,3/4 
'5,100 L 7 



(4) 



from the center, where ■ 10 7 yr is the timescale 
of the outburst, and the shock's Mach number is 
M » 2.1{AK/K 1 ) 1 ' 2 , implying that it can remain 
supersonic to a radius 



30kpc 



AK r 



(5) 



lOkeVcm 2 , 
in the power-law medium we have assumed. 

3.2. Energy-Driven Shock Heating 

After the kinetic outflow shuts off, the de- 
posited energy, E « (3 x 10 59 cvg)P^t-;, can 
drive the shock, if the front is still supersonic. 
The velocity of an energy-driven shock is v « 
/^(iJ/pr 3 ) 1 / 2 , where fp is another structure fac- 
tor analogous to fp, leading to an entropy boost 
AKp cx /}~ 5 / 3 r~ 3 . In the power-law core model, 
this entropy boost declines with radius: 

AK E (r) w 19keVcm 2 P 4 5t7/iT 5 ; i 5 / 2 

-4/3 

• (6) 



20kpc 



Making the reasonable assumption that AKp ~ 
AKp(r/rp)~ 4 / 3 , one finds that the energy-driven 
shock becomes subsonic at 



21kpc 



AK f 



lOkeVcm 2 



3/14 



T 3/14 3/7 (7 s 
J 5,100 r 7 ■ \' ) 



Notice that this form of entropy injection will 
cause an entropy inversion if the outburst is too 
powerful and too brief. That happens because an 
energy-driven shock slows down more quickly than 
an outflow-driven shock, raising the innermost gas 
to a higher entropy level than gas farther from the 
center. Boosting the entropy level in the inner 
~ 30kpc of a cluster by a uniform ~ lOkeVcm 2 
is therefore most naturally done by an outburst 
lasting several times 10 years. 

3.3. Bubble Heating 

Once the outer shock front becomes subsonic, 
the plasma bubble inflated by the AGN out- 
flow will buoyantly rise through the intracluster 
medium. Inflating the bubble has lifted some of 
the intracluster medium, increasing its gravita- 
tional potential energy by an amount correspond- 
ing to the work required to inflate the bubble. As 
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the bubble subsequently rises, the gas around the 
bubble falls. Assuming that the falling gas locally 
thermalizes that liberated gravitational potential 
energy, one finds that the bubble introduces heat 
energy amounting to 



pVbub 9 Ar = Vbub 



dP 

dr 



Ar 



(8) 



where Vbub is the bubble volume, into each ra- 
dial interval Ar (Churazov et al. 2001; Begelman 
2001). 

The amount of additional work a bubble does 
as it rises and expands, keeping PV^£ b constant, 
depends on the polytropic index 7b u b of its con- 
tents. If the bubble remains coherent as it ascends, 
it distributes the energy i?bub required to create it 
according to 



dE -Bbub 

~r — sbub 

dr 7-ini 



(9) 



where r; n j is the radius at which the bubble is 
injected and £ bu b = 2(7 bub - l)/37 bu b, equaling 
1/6 if the bubble contents are relativistic and 4/15 
if they are non-relativistic. Here we will assume 
that bubbles are injected at the radius r ss where 
the shock driven by the outburst becomes sub- 
sonic. Those rising bubbles can generate addi- 
tional entropy by causing turbulence that mixes 
low-entropy gas from the cluster core with the 
higher-entropy gas lying above it. However, we 
will not consider mixing processes here because 
they do not increase the total amount of thermal 
energy in the intraclustcr medium. 

Heating the gas by an amount Aq per parti- 
cle augments its specific entropy by A In if 3 / 2 = 
Aq/kT, implying AK = (2/3n 2 e /3 )Aq. The to- 
tal amount of bubble heating per particle from a 
single outburst is 



Aq(r) 



lim p dE 
47rr 2 p dr 



(10) 



implying that the entropy boost in the bubble- 
heated zone is Aifbub <x p~ 5 / 3 r _3 ~£ bub . Defining 
/bub = Et,uh/E to be the fraction of the input 
energy left in bubble form when the shock be- 
comes subsonic, we obtain the following expres- 
sion for bubble heating as a function of radius in 
the power-law core model 



30kpc 



"5bub 



(11) 



Like an energy-driven shock, this form of heating 
will also cause an entropy inversion if too much 
bubble energy is injected at too small a radius be- 
cause its entropy injection profile falls even more 
rapidly with radius than the profile for the energy- 
driven shock. 

3.4. Net Entropy Change 

A single AGN outburst thus generically heats 
the cluster core in three zones, an inner outflow- 
driven zone, an intermediate energy-driven zone, 
and an outer bubble-heated zone. Disturbances 
moving at the sound speed propagate to 100 kpc 
in ~ 10 8 years, during which time the inner re- 
gions relax to a new equilibrium configuration with 
slightly elevated entropy and resume cooling. In 
order to evaluate the net entropy change over a 
single cycle of heating and cooling, we will first de- 
termine the entropy losses due to radiative cooling 
during this time interval and then compare those 
losses to the entropy boost from a single outburst. 

A radiative loss rate of q = n;n e A(T)/(ni + n e ) 
carries entropy away from a plasma of ion density 
n, at the rate K = -{2/S)[m/{m + n e )]nj 3 A(T). 
Pure free-free cooling therefore removes an en- 
tropy 



AK c (r) « -9.3keVcm 2 T 54 oo 



10 kpc 



-1/6 



At 
10 8 yr 



(12) 



AXbubM 



4.0 keV Cm ^bub/bub/f; 



during a time interval At, given the simple power- 
law model we have been using for the core gas 
configuration. Notice that if outflow-driven shock 
heating can compensate for cooling at small radii, 
then heating and cooling will remain nearly bal- 
anced out to rp because of the weak dependence 
of AK C on radius. However, in the power-law core 
model, cooling will dominate heating at r > rp 
over the course of one cycle because both energy- 
driven shock heating and bubble heating fall off 
more quickly with radius than cooling. 

Thus, in order for heating to compensate for 
cooling throughout the ~ 100 kpc region where 
the cooling time is less than a Hubble time, the 
outflow-driven zone must be relatively large, im- 
plying that outbursts capable of stopping a cool- 
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ing flow indeed last several times 10 7 years. Fig- 
ure 2 illustrates this idea in the context of the 
power-law core model. Solid lines show the net 
entropy change over the course of a single cycle 
of heating and cooling for outbursts lasting 10 7 
years. Two different cases are shown, correspond- 
ing to entropy boosts of AKp = lOkeVcm 2 and 
20keVcm 2 in the outflow-driven zone. Outside 
of that, in the energy-driven zone, the boost in 
each case is Aifp(r/rp)~ 4 / 3 , where rp is given by 
equation (4) with Ts^ioo = 1. This zone extends to 
r ss as given by equation (7), and beyond that point 
bubble heating is assumed to raise the entropy by 
Aifp(r ss /rp)~ 4 / 3 (r/r ss )~ 3 / 2 , so that the entropy 
jump is continuous across r ss . The validity of that 
assumption depends on the structure factors fp, 
fs, and /bub and needs to be checked with nu- 
merical simulations. The power-law index in ra- 
dius here applies to the relativistic bubble case, 
but the non-relativistic case is virtually the same, 
with an r~ ls dependence. 

Because the time interval between heating 
episodes depends on the entropy boost at small 
radii, we choose At so that heating matches ra- 
diative losses at a radius of 2 kpc. Outside of that 
radius, the net entropy gain is several keV cm 2 out 
to about rp, where there is a sharp turnaround. 
Beyond about 30 kpc the net entropy change 
is negative, even for a central entropy boost of 
20 keV cm 2 . The dashed lines show how the pic- 
ture changes when an outburst lasts long enough 
for rp to equal r ss as given by equation (5). Then 
the outflow-driven zone extends much farther, and 
in the case of a 20 keV cm 2 boost, heating can 
exceed cooling out to ~ 100 kpc. The time re- 
quired for an outflow-driven shock to reach r ss is 
2.4 x 10 7 T 5 " 1 1 qq yr in the 10 keV cm 2 case and 
4.3 x 10 7 T 5 ; 1 1 / O 2 yr in the 20 keV cm 2 case. 

A slightly more refined picture emerges when 
we take the core structure at the time of the out- 
burst to be that of the pure cooling model of Voit 
et al. (2002), corresponding to the solid line in 
Figure 1. In that model of the core, gas density 
starts to drop with radius more rapidly than r~ 4 
outside of ~ 30 kpc, gradually increasing the ef- 
ficiency with which all forms of heating increase 
the entropy. Figure 3 shows models in which the 
entropy boost from the outburst is a continuous 
piecewise power law that is oc (rp) -4 / 3 in the 
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Fig. 2. — Net entropy change AK nc t during a single cycle 
of heating and cooling given the power-law model of core 
structure. The upper set of lines (blue) shows models in 
which each AGN outburst heats the intracluster medium 
at a radius of 2 kpc by 20 keV cm 2 . The lower set of lines 
(red) shows models in which each AGN outburst heats the 
intracluster medium at a radius of 2 kpc by 10 keV cm 2 . 
The cooling time after each outburst is normalized to give 
zero net entropy change at that radius. Solid lines show 
the net entropy change for an outburst lasting 10 7 years 
that deposits heat according to the piecewise power law 
described in the text. Dashed lines show the net entropy 
change for outbursts lasting long enough for the heating to 
remain outflow-driven until the shock front becomes sub- 
sonic. The maximum net entropy gain occurs at the radius 
rp the shock has reached at the time the outburst ceases. 
Only the longer outbursts are able to prevent net cooling at 
> 30 kpc, suggesting that each AGN outburst must remain 
at ~ 10 45 erg s — 1 for several times 10 7 years in order to shut 
off the cooling flow while maintaining a quasi-steady inner 
entropy profile. 

outflow-driven zone, oc p~ 5 / 3 r~ 3 in the energy- 
driven zone, and oc p~ 5 /3 r -3-£ b ub m ^be bubble- 
heated zone. Boundaries between the zones are 
the same as those used for Figure 2. The entropy 
boost is normalized so that AKp at 2 kpc takes 
values of 10, 15, 20, and 30 keV cm 2 , and cooling 
is assumed to match heating at that radius. Solid 
lines show models for outbursts lasting 10 7 years, 
and dashed lines show models for outbursts long 
enough for the outflow-driven zone to reach the 
radius at which the shock becomes subsonic. 

The simple models illustrated in Figure 3 
suggest that episodic outbursts that impart ~ 
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Fig. 3. — Net entropy change AK nB t during a single cy- 
cle of heating and cooling when the pure cooling model of 
Voit et al. (2002) is used for the pre-outburst structure 
of the core. Each set of lines shows two models with the 
same value of entropy input at 2 kpc. From top to bot- 
tom, these values are 30 keV cm 2 (magenta), 20 keV cm 2 
(blue), 15 keV cm 2 (green), and 10 keV cm 2 (red). Cooling 
after each outburst is normalized to give zero net entropy 
change at 2 kpc. Solid lines show the net entropy change 
for an outburst lasting 10 7 years. Dashed lines show the 
net change for outbursts lasting long enough for the heat- 
ing to remain outflow-driven until the shock front becomes 
subsonic. Interestingly, an outburst that raises the cen- 
tral entropy level to 15 keV cm 2 marginally compensates 
for cooling at all radii over the course of a cycle, suggest- 
ing that the characteristic minimum central entropy level 
of ~ lOkeVcm 2 is determined by the minimum outburst 
required to halt a cooling flow. 



15keVcm 2 to the central regions of a cluster and 
last until the shock reaches ~ 30 kpc can com- 
pensate for radiative cooling throughout an entire 
cluster. In all the models, heating balances cool- 
ing from 1 kpc to 5 kpc with a high precision that 
is likely to be coincidental, given the crudity of the 
model. Outside of 5 kpc, the net entropy gain per 
cycle rises out to the outer radius of the outflow- 
driven zone, beyond which the net entropy gain de- 
clines because the shock decelerates more quickly 
there. The sharp cusp appearing at the transi- 
tion to the energy-driven zone is an artifact of the 
piecewise power-law model for heating; a more re- 
alistic model would have a more gradual turnover. 
Beyond ~ 50 — 100 kpc, in the bubble-heated zone, 
the net entropy change climbs again with radius 



because the gas density profile gradually steepens 
there, approaching p oc r~ 2 beyond the cluster 
core. The entropy gain outside the core from bub- 
ble heating is therefore Aifbub oc r 1 ^ 3 ^^^ , which 
is nearly constant with radius. However, bubble 
heating operates out there only if the bubbles can 
remain coherent over many pressure scale heights, 
requiring some sort of stabilizing mechanism to 
ward off hydrodynamic shredding (e.g., Bruggen 
& Kaiser 2001). 

One intriguing possibility suggested by Figure 3 
is that the characteristic minimum entropy level of 
~ lOkeVcm 2 in cluster cores is determined by the 
minimum AGN outburst required to halt a cooling 
flow. All the cluster entropy profiles in Figure 1 
flatten at roughly this level even though the clus- 
ters range over a factor more than three in temper- 
ature. Our analysis indicates that AGN outbursts 
unable to heat the central regions to this level are 
also not powerful enough to push outflow-driven 
shocks to ~ 30 kpc, the necessary inner radius for 
either energy-driven shocks or bubble heating to 
compensate for radiative cooling in the outer parts 
of a cluster's core (~ 50 — 100 kpc). This feature of 
our simple model suggests that the AGN feedback 
mechanism may have become tuned over time so 
that cooling is kept marginal (see Omma & Binney 
2004, for an example of such a mechanism). 

An alternative possibility is that the consis- 
tency among the inner entropy levels in the sample 
of Donahue et al. (2005) is a selection effect, be- 
cause clusters with the lowest central entropy lev- 
els will appear to have the strongest cooling flows. 
It could be that outbursts of kinetic energy from 
a cluster's central AGN outbursts vary greatly in 
power, sometimes imparting entropy boosts sev- 
eral times greater than lOkeVcm 2 . In that case, 
however, one would expect a large fraction of clus- 
ters to have elevated central entropy levels because 
it would take much longer for radiative cooling to 
reduce the central entropy in those clusters. 

3.5. Multiple Outbursts 

Now we will briefly consider the evolution of 
the intracluster medium over time in the context 
of this simple framework for AGN heating. In gen- 
eral, one would not expect precise balance between 
AGN heating and radiative cooling throughout the 
cluster core because the conditions governing ac- 
cretion at < 1 kpc are not closely linked to those 
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at ~ 100 kpc. A common feature of the out- 
burst models presented here is that heating ex- 
ceeds cooling in the outer parts of the outflow- 
driven zone, and that excess can persist well be- 
yond the outflow-driven zone if bubble heating is 
effective. One therefore expects entropy to accu- 
mulate over time in the ~ 30 — 100 kpc range, 
if the outbursts are strong enough to prevent net 
cooling everywhere outside of ~ 1 kpc. 

Evidence for such an accumulation of entropy 
is present in Figure 1. The entropy profiles shown 
there tend to exceed the pure-cooling model by 
~ 50keVcm 2 at ~ 100 kpc. The models we have 
constructed here are too crude for a detailed com- 
parison with that excess to be fruitful. Statistical 
variations in the duration and kinetic power of the 
outbursts will cause the size of the outflow-driven 
zone to differ from outburst to outburst, and the 
shape of that zone is likely to be elongated, ow- 
ing to the asymmetry of the outbursts. Neither of 
those features are represented in our spherically- 
symmetric, single-outburst models. 

Nevertheless, a generic prediction of this frame- 
work is that highly energetic outbursts should 
limit the total accumulation of entropy at large 
radii over the course of several billion years. From 
equation (1) for the cooling time from a given en- 
tropy level, one infers that the maximum number 
of outbursts that can occur over a time period t is 

m / AK P \ - 3/2 ( kT \ 

x (fife) • (13) 

if each outburst imparts an entropy boost AKp 
in the outflow-driven zone. Thus, several tens 
of ~ lOkeVcm 2 outbursts can occur, but only a 
few at the ~ 50keVcm 2 level. If the excess en- 
tropy introduced at large radii by a single out- 
burst is oc AKp, then the accumulated excess 
entropy arising from a series of such outbursts 

— 1/2 

is cx AK p ' . However, an adequate calibration 
of the relation between the central entropy input 
from an outburst and the corresponding excess en- 
tropy it produces will require high-resolution nu- 
merical simulations. 

In principle, it may be possible to use the ex- 
cess entropy observed in the outer regions of clus- 
ters to constrain the outburst history. Figure 3 in- 



dicates that each outburst incrementally deposits 
several keV cm 2 at ~1 Mpc, where cooling is neg- 
ligible, as long as bubble heating still operates 
at such large radii. For example, suppose a se- 
ries of ~ 15keVcm 2 outbursts were to introduce 
<~ 5keVcm 2 at <~ 1 Mpc with each outburst, as 
suggested by Figure 3. Roughly thirty such out- 
bursts could occur in a ~ 5 keV cluster over a pe- 
riod of 5 Gyr, according to equation (13), leading 
to an accumulated excess entropy ~ 150keVcm 2 
at large radii. The entropy excesses seen in the 
outer parts of < 5keV clusters (Ponman ct al. 
2003) might therefore result from the cumulative 
effects of bubble heating during the cluster's life- 
time (see also Roychowdhury et al. 2004). 

In even smaller systems with temperatures < 
1 keV, excess entropy that accumulates at > 
100 kpc may come to dominate the gravitation- 
ally produced entropy. Because these systems 
have smaller core radii than clusters, AGN out- 
bursts can penetrate more easily to radii beyond 
~ 30 kpc, depositing more entropy there. Thus, 
the "entropy deficit" observed in the outskirts of 
groups by (Mahdavi et al. 2005) might not be 
an actual deficit in the gravitationally produced 
entropy but rather just a deficit with respect to 
a large entropy excess of non-gravitational origin 
that has built up near the core radius. 

3.6. Outbursts of Sustained Power 

The dramatic heating event observed in MS0735 
+7421 by (McNamara et al. 2005) inspires us 
to add one more wrinkle to this framework for 
AGN heating. So far, we have considered only 
outbursts of ~ 10 45 ergs _1 lasting significantly 
less than 10 8 years, but MS0735+7421 appears 
to have been sustaining a kinetic power output 
~ 1.7 x 10 46 ergs _1 for ~ 10 8 years. The observed 
central entropy level of ~ 30 keV cm 2 is consistent 
with a constant-power outburst at this level, given 
a structure factor fp ~ 0.5, and corresponds to a 
central cooling time ~ 1 Gyr at the central tem- 
perature of 3 keV. In those respects MS0735+7421 
is similar to the two radio-quiet cooling-core clus- 
ters observed by Donahue et al. (2005b), which 
also have ~ 30 — 50keVcm 2 entropy levels and 
~ 1 Gyr central cooling times, suggesting that 
these clusters are what MS0735+7421 will look 
like a few hundred million years from now. If 
that is the case, then instances of powerful, sus- 
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taincd outflows from AGNs in clusters might not 
be all that rare and should be incorporated into 
the overall framework for AGN heating. 

How much longer the energetic outflow in 
MS0735+7421 will persist is anyone's guess, but 
its effects on the cluster's overall entropy profile 
are likely to be profound and lasting. The outflow 
has already propagated to a radius ~ 300 kpc, 
leaving the cluster core and entering the region 
where gas density is steepening to a p oc r~ 2 den- 
sity profile. An interesting feature of an outflow- 
driven shock in a medium with that density distri- 
bution is that it propagates at constant velocity, 
generating an entropy boost that is proportional 
to the original entropy, i.e. AK/K is constant. 
A sustained, high-power outflow therefore raises 
the outer entropy profile of a cluster by a constant 
factor, preserving its shape. 

That sort of entropy boost is particularly 
interesting in light of XMM-Newton observa- 
tions showing that the outer entropy profiles of 
low-temperature clusters have shapes similar to 
those of hotter clusters but have unexpectedly 
large normalizations, indicating that some non- 
gravitational process has elevated the entropy 
by a constant factor throughout the intracluster 
medium (Pratt & Arnaud 2003, 2005). One sug- 
gested possibility for causing this global entropy 
boost is that early galactic winds smoothed the in- 
tergalactic medium of the protocluster before the 
cluster formed, thereby amplifying the amount of 
entropy generated in the cluster's accretion shocks 
(Voit et al. 2003; Ponman et al. 2003). However, 
simulations of entropy amplification suggest that 
supernova-driven winds have trouble generating 
enough smoothing to account for the observed en- 
tropy excesses (Borgani et al. 2005). Smoothing 
by AGNs might be required, or maybe it is pos- 
sible for the central AGN to produce the entire 
entropy excess through large outbursts alone. 

In order to assess the large-scale impact of heat- 
ing events like that in MS0735+7421, we can con- 
sider an idealized cluster model in which the clus- 
ter atmosphere has p oc r~ 2 everywhere. Taking 
the baryon fraction to equal the global value gives 
pr~ 2 w 3.2 x lO^gcm- 1 (7cT/5keV), where T is 
now the mean temperature of the cluster. The en- 
tropy boost produced in that density profile by a 



sustained outburst of kinetic power output is 

AK M m p /g, /P kin \ 2/3 

K 3 • 4 2 / 3 fcT ^ pr 2 ) 

where = Pkin/10 46 ergs -1 . In other words, 
the outburst event in MS0735+7421 is destined 
to raise its global entropy profile by > 10%, de- 
pending on the structure factor fp. Further- 
more, lower-temperature clusters are subject to 
even larger global boosts, implying that strong, 
sustained outbursts of AGN heating could be a 
substantial contributor, maybe even the dominant 
contributor, to the entropy excesses observed at 
large radii in low-temperature clusters. 

4. Summary 

In this paper we have outlined a framework for 
AGN heating of cluster cores inspired by the ob- 
served core entropy profiles of cooling- flow clusters 
(Donahue et al. 2005a). We have focused primar- 
ily on entropy rather than on density or temper- 
ature, because entropy is the quantity most di- 
rectly affected by heat input and radiative losses. 
Tracking entropy changes allows one to predict, at 
least qualitatively, how an outburst of AGN en- 
ergy becomes distributed with radius in the in- 
tracluster medium. Applying these simple pre- 
dictions to the observed distributions of entropy, 
temperature, and density in cluster cores indicates 
how the entropy at each radius should change over 
a single cycle of outburst heating and radiative 
cooling. The picture we arrive at from consid- 
ering the time-dependent behavior of these core 
entropy profiles has much in common with other 
AGN-heating models that were inspired by other 
aspects of clusters with cool cores, especially the 
X-ray cavities often found there. Here we sum- 
marize the AGN-heating framework that emerges 
from considerations of intracluster entropy. 

The key features of the core entropy profiles in 
cooling-flow clusters, described in Section 2, are as 
follows. 

1. The profiles observed by (Donahue et al. 
2005a) are very similar to one another when 
plotted as functions of radius in unsealed 
physical units, even though the clusters 
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range over a more than a factor of three in 
temperature. Power-law fits to these profiles 
show that K oc r a with a « 0.9 — 1.3, and 
the normalizations of the profiles are consis- 
tent with earlier indications that entropy at 
a fixed fraction of the virial radius is oc X^ 3 , 
where is the characteristic temperature of 
the cluster's halo (Ponman et al. 2003; Voit 
& Ponman 2003; Pratt & Arnaud 2003; Pif- 
faretti et al. 2004) . A power-law model with 
K{r) « 150keVcm 2 (r/lOOkpc), T oc r 1 / 3 , 
and p oc r _1 is a reasonable approximation 
to the gas configuration from ~ 10 kpc to 
- 100 kpc. 

2. In the outer regions of these cluster cores, 
the entropy profiles appear to approach the 
pure-cooling cluster model of Voit et al. 
(2002), in which a cluster whose initial 
configuration is determined by hierarchi- 
cal structure formation is allowed to cool 
for a Hubble time without feedback. The 
pure-cooling model is apparently a bounding 
lower envelope to the allowed set of entropy 
profiles. 

3. The inner regions of these profiles flatten out 
at an entropy level ~ lOkeVcm 2 . Further- 
more, simply adding an entropy pedestal of 
lOkeVcm 2 to the pure-cooling model at all 
radii reproduces the observed shape of these 
profiles at radii from ~ 1 kpc to > 30 kpc. 

4. This inner entropy minimum corresponds 
to a cooling time ~ 10 8 years, suggest- 
ing that the heating mechanism producing 
the entropy pedestal is episodic on roughly 
this timescale (see also David et al. 2001; 
Kaiser & Binney 2003). If AGN heating is 
indeed episodic, then the relatively strong 
power-law dependence of cooling time on 
entropy implies that the central entropy 
level remains near its maximum value during 
most of the time period between feedback 
episodes. Thus, typical AGN heating events 
probably impart ~ 10 — 15keVcm 2 in the 
central parts of the core. 

5. Because the entropy pedestal introduced 
by a heating episode persists over most of 
the cycle, the central temperature gradi- 
ent should remain approximately flat within 



~ 10 kpc, which corresponds to a tempera- 
ture that is approximately one-half the value 
at ~ 100 kpc, given the T oc r 1 / 3 temper- 
ature gradient outside of 10 kpc (see also 
Kaiser & Binney 2003). In this configura- 
tion, the lower bound to the gas tempera- 
tures observed in cooling-flow clusters arises 
in a quasi-steady hydrostatic configuration 
with an inner entropy floor and reflects the 
potential-well depth in the region where the 
entropy distribution flattens. 

In addition to these features of the entropy pro- 
files, the strong iron abundance gradients seen in 
the cores of cooling-flow clusters provide another 
important constraint on the heating mechanism: 
it cannot induce too much mixing in the core. 
Otherwise, it would wipe out the observed gra- 
dients (e.g., Bohringer et al. 2002). 

Taken as a whole, these features of the entropy 
structure in cluster cores point toward a specific 
framework for AGN heating, independent of ob- 
servations of X-ray cavities (see §3). 

1. An AGN outburst with a constant kinetic 
power output naturally produces a con- 
stant entropy pedestal if it propagates into 
a medium with p oc r _1 . Raising the inner 
entropy by ~ lOkeVcm 2 requires outbursts 
of ~ 10 45 ergs _1 , given the values of rp(r) 
observed in the cores of cooling-flow clusters. 
This mode of heating is particularly attrac- 
tive because a boost in entropy that is con- 
stant with radius does not produce entropy 
inversions, which are not generally observed 
and would induce convection, mixing, and 
reduction of the central iron gradient. 

2. In order to produce an entropy pedestal ex- 
tending to ~ 30 kpc, as suggested by the 
observed profiles, an outburst sufficient to 
raise the central entropy by ~ lOkeVcm 2 
needs to last several times 10 7 years. Shorter 
bursts pump too much outburst energy into 
too small a volume, leading to entropy in- 
versions, because both energy-driven shocks 
and bubble heating produce entropy boosts 
that decline with radius in a medium where 
p oc r _1 . 

3. The leading shock of an outburst that raises 
the central entropy by ~ lOkeVcm 2 and 
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lasts for ~ 3 x 10 7 years will become sub- 
sonic and decay to sound waves in the neigh- 
borhood of ~ 30 kpc from the center. At 
that point the hot bubble of ejected material 
that was driving the shock can become buoy- 
ant. If the intraclustcr medium efficiently 
thermalizes the work done by that bubble 
as it rises, then it can offset cooling at large 
radii in the core (Churazov et al. 2000, 2001, 
2002). Small bubbles can be just as effective 
as large ones, as long as the total amount of 
initial bubble energy is the same (Begelman 
2001; Ruszkowski & Begelman 2002). One 
can therefore compute the entropy injected 
at all radii by an outburst, given the size of 
the central entropy boost and the duration 
of the outburst, under the assumption that 
the bubbles remain coherent as they rise. 

4. Assuming that the outbursts are regulated 
by the need for heating to balance radiative 
cooling at small radii over the course of an 
outburst cycle determines the net entropy 
change at all radii from a single cycle. The 
size of the central entropy boost determines 
the number of outburst cycles that can oc- 
cur during the life of a cluster, and therefore 
regulates the accumulation of excess entropy 
at larger radii. 

5. In order for bubble heating to exceed cool- 
ing in the outer parts of the core over the 
course of a cycle, the inner radius of the 
bubble-heated region must be > 30 kpc. 
That requirement places a lower limit ~ 
10 — 15keVcm 2 on the size of the central 
entropy boosts. Weaker outflows do not ex- 
tend the base of the bubble-heated zone to 
a sufficiently large radius. This feature of 
cluster cores may explain the ubiquity of the 
<~ lOkeVcm 2 pedestal. Perhaps the feed- 
back mechanism has become tuned over time 
to produce outbursts with power output and 
duration sufficient to ward off cooling catas- 
trophes. 

6. Outside the cluster core, the gas density pro- 
file steepens to p oc r~ 2 , thereby increasing 
the efficiency of entropy production. Sus- 
tained outbursts with constant kinetic power 
output can therefore produce entropy boosts 
for which AK/K is constant with radius, 



raising the previous entropy profile by a con- 
stant factor without changing its shape. An 
outburst power of ~ 10 46 ergs _1 , such as 
that observed in MS0735+7421, can raise 
the outer entropy profile of a 5 keV cluster by 
> 10%, and the boosts in lower-temperature 
clusters can be substantially greater. Strong 
AGN outbursts can therefore potentially ac- 
count for the elevated outer entropy profiles 
of low-temperature clusters (see also Roy- 
chowdhury et al. 2004). 

Interestingly, the picture that emerges from 
consideration of the core entropy profiles alone 
calls for outbursts that deposit a total energy 
~ 10 60 erg every ~ 10 s years, corresponding to a 
time-averaged AGN power output of several times 
lO^ergs" 1 with a duty cycle - 20 - 30%. Di- 
viding this energy between two bubbles, each re- 
quiring an energy of 4PV to produce, implies that 
PV ~ 10 59 erg in each bubble. All of these charac- 
teristics are remarkably consistent with the prop- 
erties of AGN outbursts inferred from X-ray cavi- 
ties (e.g., Birzan et al. 2004) 

Further progress within this framework will re- 
quire numerical simulations to determine the ac- 
tual entropy boost imparted by a single outburst 
as a function of radius, given an initial config- 
uration for the core. There has already been a 
substantial amount of work in this area, ranging 
over a wide volume of parameter space in kinetic 
power, outburst duration, and core configuration 
at the time of the outburst (recent papers include 
Dalla Vecchia et al. 2004; Omma et al. 2004; Zanni 
ct al. 2005). The analysis we have presented in- 
dicates the most fruitful area of this parameter 
space to explore in future work, while strength- 
ening the case for episodic AGN outbursts as the 
mechanism that limits cooling and condensation 
in cooling-flow clusters at the present time and 
therefore the growth of the largest galaxies in the 
universe. A complete picture of the intracluster 
medium will also need to account for earlier pro- 
cessing of a cluster's intergalactic gas by cooling 
and supernova feedback during the epoch of galaxy 
formation, whose effects can be amplified by accre- 
tion shocks as the cluster forms (Voit et al. 2003; 
Borgani et al. 2005). 
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